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The equilibrium shapes of InAs quantum dots (i.e., dislocation-free, strained islands with sizes 
> 10, 000 atoms) grown on a GaAs (001) substrate are studied using a hybrid approach which 
combines density functional theory (DFT) calculations of microscopic parameters, surface energies, 
and surface stresses with elasticity theory for the long-range strain fields and strain relaxations. 
In particular we report DFT calculations of the surface stresses and analyze the influence of the 
strain on the surface energies of the various facets of the quantum dot. The surface stresses have 
been neglected in previous studies. Furthermore, the influence of edge energies on the island shapes 
is briefly discussed. From the knowledge of the equilibrium shape of these islands, we address 
the question whether experimentally observed quantum dots correspond to thermal equilibrium 
structures or if they are a result of the growth kinetics. 



I. INTRODUCTION 

In recent years, the study of growth conditions and 
electronic properties of quantum dots has attracted 
nificant attention in basic science and technology.E 
Quantum dots are small three-dimensional islands of a 
low-band-gap semiconductor (e.g. In^jGai-a^As) which 
are enclosed in a wide-band-gap semiconductor matrix 
(e.g. GaAs). Provided the bands of these islands and of 
the host are appropriately aligned, the valence and con- 
duction bands produce a confinement potential for the 
holes in the valence band and the electrons in the con- 
duction band. If these islands are small enough, they will 
behave like big artificial atoms with discrete energy lev- 
els. Thus, the recombination spectrum of a single quan- 
tum dot consists of a single sharp line with a practically 
not measurable temperature broadening. Quantum dots 
may be used for new types of devices, as for example 
a single-electron transistor or cellular automata. Fur- 
ther examples are semiconductor lasers, where the wave 
length of the emitted light is determined by alloy com- 
position and the size and shape of the dots. Indeed, such 
lasers have been build in the laboratory.BEl The required 
size of the quantum dot is dictated by the condition that 
the energy separation of the quantized electronic levels of 
the dots should be about 0.1 — 0.2 eV, so that they are not 
populated at room temperature. And the dot has to be 
sufficiently large that at least one bound level exists; for 
too small islands this is not the case.B For GalnAs these 
conditions imply that the width of the island is between 
50 and 200 A (or 20 — 80 atoms), which means that we 
are dealing with objects built of 1,000 — 60,000 atoms. 
For lasers it is necessary to have many dots (~ lO^^cm"^) 
and these should all have nearly the same size and shape. 



so that all dots emit light at practically the same wave 
length; in a quantum dot laser the width of the line is 
determined by the size and shape fluctuations of the en- 
semble of quantum dots which implies that the size and 
shape uniformity of quantum dots is critical to these ap- 
plications. 

Already in 1990 it had been observed that dislocation- 
free, strained islapds form by itself when InAs is de- 
posited on GaAs,EJ and since 1994 several groupg3~u have 
shown that a range of growth parameters exists at which 
quantum dots assemble themselves with the desired and 
tunable size and a rather narrow size distribution. The 
mechanism giving rise to this self-assembly of the dots is 
still not understood. 

Often the formation of quantum dots is explained in 
terms of a thermal-equilibrium picture where the system 
assumes the state of lowest free energy: Islands form, 
instead of a strained, epitaxial film, because the gain 
of elastic relaxation energy (possible in an island) over- 
compensates the cost due to the increased surface energy 
(a three-dimensional island has a larger surface than a 
two-dimensional film). Typically such island do not form 
immediately on the substrate but on top of a wetting 
layer (see for example Ref. ||). For InAs quantum dots 
on GaAs (001) this wetting layer has a thickness of about 
1.5 monolayers. When this thermal equilibrium picture 
applies the growth mode giving rise tp islands is called 
the Stranski-Krastanov growth modeB 

OtheiL^ authors have emphasized the role of kinetic 
effects.E£rEl Dobbs et alLn studied the formation of is- 
lands using self-consistent rate equations. Their rate the- 
ory is designed to predict reliably average quantities and 
they found that their island density in dependence of the 
coverage is in good agreement with experiment. But they 



could not calculate the island size distribution explicitly. 

Becausepp|r^veral experimental as well as theoreti- 
cal papersQMIlil it is assumed that thermal-equilibrium 
theory is applicable to describe and understand the for- 
mation of quantum dots, we performed calculations of 
the equilibrium structure to find out whether or not this 
agrees with results of growth experiments. The calcu- 
lations are done for dislocation-free InAs islands epitax- 
ially grown on GaAs (001). Knowing the equilibrium 
shapes is indeed important, because under certain con- 
ditions thermal equilibrium will be reached. And in gen- 
eral, if the quantum-dot shape observed in experiments 
deviates from the equilibrium shape, one has to conclude 
that equilibrium thermodynamics is not adequate to de- 
scribe the island formation and the assumed shape and 
size distribution. We note in passing that the growth of 
experimental quantum dots is likely affected^lso by (un- 
wanted) alloying of the dot and the matrix)l3 by entropy 
effects, and for high concentrations of quantum dots also 
the island-island interaction has to be considered. These 
effects are neglected in the present calculations. 

The InAs quantum dots of interest contain more than 
1,000 atoms. We consider it unwise to evaluate the total 
energy of such a system by a direct density functional cal- 
culation. Therefore a hybrid method was developed that 
provides the same result as a direct approach, but at 
much less computational costs. Even more importantly, 
this method exhibits with greater clarity the underlying 
physical mechanisms. In brief, the approach is summa- 
rized as follows. The total energy of a large, isolated 
quantum dot is given by 
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The leading terms are the elastic relaxation energy and 
the sum over the surface energies of the surface facets. 
Both quantities depend sensitively on the quantum-dot 
shape. The surface reconstructions, the surface energies 
and their strain dependence are calculated by DFT and 
analyzed as a function of the atomic chemical potential. 
If the size of an island is bigger than about 1,000 atoms 
it turns out that the strain fields and elastic energies 
are well described byj-elasticity theory because they fol- 
low the scaling laws.ll3 We therefore evaluate the long- 
range strain relaxation in the quantum dot and in the un- 
derlying substrate by elasticity theory applying a finite- 
element approach. The approach permits the systematic 
investigation of almost any island shape. Fig. |l| displays 
the importance of the different energy contributions and 
their scaling with island size. These results refer to an 
isolated, pyramidal shaped quantum dot, but essentially 
the same behavior is found for other island shapes. The 
elastic energy of a relaxed quantum dot, compared to 
the energy of a two-dimensional epitaxial film scales lin- 
early with the volume (or the number of atoms) of the 
quantum dot. This analytical scaling holds true for the 
elastic energy as long as the islands contain more than 
1,000 atoms. The single facets of the island have to be 
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FIG. 1. Energy contributions of a single, pyramidal 
shaped, dislocation free island compared to the energies of 
an epitaxial (strained) film. The elastic energy relief E^^^^^^^ 
due to the partial strain relaxation, the surface-energy con- 
tribution E^^^^^^^ due to the increased surface area, and the 
contribution due to the various edges E^^^^ are shown. The 
total energy displays the sum of all these contributions. The 
applied approach (see text) is valid for island sizes bigger than 
1,000 atoms. The elastic constants used for the strain relax- 
ation are given in Tab. Q The surface energies stresses are 
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those of Tab. 
meV/A. 



^. For the edge energy we used 7^^se 



larger than ~ 16 atoms, so that the reconstructions on 
the facets are not suppressed and therefore the surface 
energy scales then with the area. This implies it scales 
with the volume to the power 2/3. Surface energy is 
a cost, and therefore the contribution is positive. Also 
shown is the contribution from the edges which also is a 
cost, though a rather small one. Obviously, the sum over 
the energies of the edges scales with the volume to the 
power 1/3. So that the edge energies provide the scaling 
relation they should be larger that ~ 4 atoms. The main 
uncertainty in using the above approach is that there 
are no edge energies known until now and they only can 
be estimated. Furthermore, for small islands the atomic 
structure on the side facets might not reconstruct. 

This method of calculating the total energy of 
dislocation-free, strained, jielaxed islands (as described 
by Eq. |^) was used beforeliZl'ES but so far some approx- 
imations were implied which will now be dropped. The 
approximations were: i) the elastic properties of the InAs 
islands were taken identical to those of the GaAs sub- 
strate, ii) the influence of the surface strain on the sur- 
face energies of the various facets was neglected, in) the 
influence of the edge energies was not discussed. The im- 
proved treatment reported in this paper required elabo- 
rate calculations (this applies in particular to the surface 
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stresses), but our results fully confirm the earlier con- 
clusions. The obtained quantitative differences to the 
previous work are small. 

In the following Section we discuss some results of the 
finite element calculations. Then, Section III presents 
the calculations of surface energies and surfaces stresses 
for the low-index surfaces: (001), (110), (111), and (111). 
For the (110) surface we also give the results for the first 
derivative of the surfaces stress with respect to strain. 
Combining the results for elastic and surface energies we 
obtain the total energy of the islands, and from the total 
energies of all different island shapes we derive the ener- 
getically favorable island shape for a given volume. This 
analysis is done in Section IV, where we also discuss the 
influence of the edge energies. 



II. LONG-RANGE STRAIN RELAXATION IN 
THE ISLAND AND THE SUBSTRATE 

We have calculated the elastic energy within the con- 
tinuum theory. The experimental elastic moduli (see 
Tab. 1) are used to describe the elastic properties of both 
the substrate and the island. A finite element approach 
is applied to solve the elasticity problem. Both the is- 
land and a sufficiently thick slab (240 A for a quantum 
dot volume of 2.88 x 10^ A^) representing the substrate 
are divided into small irregularly shaped tetrahedra. The 
displacement field is tabulated on the vertices of this par- 
titioning. Within each tetrahedron the linear interpola- 
tion of the displacement field is uniquely determined by 
the values at the four corners of the tetrahedron. The 
total elastic energy is calculated by summing the elastic 
energy density within each tetrahedron, which is a func- 
tion of local strain, times the volume of the unstrained 
tetrahedron over all tetrahedra. This expression is itera- 
tively minimized with respect to the displacement field. 
This procedure is repeated for finer and finer partitioning 
of the volume and the results are finally extrapolated to 
fineness equal to zero. 

To obtain the elastic energy of truncated islands we use 
a simple analytic approximation to avoid repeating the 
full finite-element calculation. This analytic expression 
is based on the scaling law for the elastic energy and on 
the fact that the tops of the pyramidal islands are almost 
completely relaxed. The islands relaxes about 50 %, if 
the substrate is kept fixed. Additional 15 % are gained if 
both island and substrate are relaxed. In Fig. || the strain 
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TABLE I. The e^^perimental elastic moduli cn, ci2 and C44 
of GaAs and InAsH 



FIG. 2. The trace of the strain tensor for (a) a pyramidal 
and (b) a truncated island on the (010) cross sections through 
the islands. Note that the actual calculation has been carried 
out for a much thicker slab. 



field of a pyramidal and a truncated island are compared. 
For both shapes the trace of the strain tensor is shown 
on a (010) cross section plane. The top of the pyramidal 
island is almost fully relaxed and therefore the elastic 
energy is almost completely stored in the base. From 
this observation and the scaling law of elastic energy with 
volume, one easily can derive the analytic approximation 
for the elastic relaxation energies of truncatGd, pyramids 
which was proven to be sufficiently accurate.ll3 



III. SURFACE ENERGIES, SURFACE STRESSES, 
AND THEIR FIRST DERIVATIVES 

The InAs surface energies and surfaces stresses are 
calculatedE2l using density-functional theory and the 
local-density approximation for the exchange-correlation 
energy functional. Ell We use ab ^^^Pri norm-conserving, 
fully separable pseudopotentials.t3~Ej The wave func- 
tions are expanded into plane waves with an energy cutoff 
of 10 Ry. The k-summation is done by using a uniform 
Monkhorst-Pack meshE^ with a density equivalent to 64 
k-points in the whole (1x1) surface Brillouin zone of 
the (100) surface. To obtain the absolute surface energies 
for (111) and (111) orientations we employ the energy 
density formalism introduced by Chetty and Martin.EEl 
Corresponding calculations were done before for GaAs 
and are described in Ref. ^ As the In A^ surface re- 
constructions are similar to those of GaAs,E3 we choose 
the same candidates for the low-energy surface struc- 
tures. Indeed, we find the same behavior, except that 
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TABLE 11. Surface energies 7 and surface stresses ax, Cy 
for InAs surface reconstructions in equilibrium with bulk As. 
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FIG. 3. Atomic structure models for the different InAs sur- 
faces, top and side views. Filled and open circles denote In 
and As atoms, respectively. 
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FIG. 4. InAs surface energies of the (110), (100), (111), 
and (111) surface orientation as a function of the As chemical 
potential. The thick lines highlight the calculated surface 
energies of the reconstructions of lowest total energy. 



As-rich reconstructions like the As terminated (110) sur- 
face are energetically unfavorable and not thermodynam- 
ically stable. 

The relaxed atomic surface geometries of the equilib- 
rium structures are displayed in Fig. ||. In Fig. ^ the 
surface energies are shown as function of the As chemical 
potential. The left and right vertical dashed lines de- 
note a In and As-rich environment, respectively. For the 
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(100) orientation the a(2 x 4) and (52(2 x 4) reconstruc- 
tions have the lowest surface energy. Both the (110) and 
(111) surface energies are independent of the As chemi- 
cal potential. The relaxed (1x1) cleavage surface is the 
stable reconstruction for the (110) orientation and the In 
vacancy structure for the (111) orientation. On the (111) 
surface an As trimer reconstruction forms in As-rich envi- 
ronment. In In rich environment the In adatom structure 
is energetically preferred. We note that a (V^ x ^/T9) 
structure had been observedin the case of GaAs by scan- 
ning tunneling microscopy .Ej However, we have not yet 
carried out calculations for this reconstruction, which 
would be rather expensive due to the large unit cell. 

Because epitaxial growth is often performed under As- 
rich conditions, we present in Tab. || the surface energies 
for /iAs = MAs(buik)7 ^'^'1 MAs IS taken at the value of the 
right dashed line of Fig. ^. Furthermore, we have cal- 
culated the surface stress for the reconstructions stable 
under As-rich conditions, and we calculated the strain 
derivatives of the stress. These results are required in 
order to obtain the corrections of the surface energy for 
strained systems. The surface energy of a strained sur- 
face defined with respect to the area of the undeformed 
surface is given by 



^, strained ^, , \ ^ 

7 =7 + 



, (2) 



where 7 is the unstrained surface energy, (Jij the surface 
stress tensor, Cij the strain tensor, and Sijki the tensor 
of second order stresses. The calculation of the first or- 
der surface stress is done as follows. We calculate the 
surface energy of a slab for various lattice constants in 
the range of ±4 %. The strained surfaces do not have to 
be relaxed again after straining because the relaxation 
energy is of second order in the strain. The energies 
of the strained surfaces are fitted to a polynomial from 
which we extract the linear coefficient of surface energy 
as a function of the strain. In all calculations we find 
that the components of the surface stress tensor are ten- 
sile. Compared to the Sij-QOO) surface its value has the 
same order of magnitudeBj For the (110) surface ax and 
ay denote components parallel to [001] and [110] respec- 
tively. For the (111) and (111) surface the surface stress 
tensors are isotropic due to the three fold symmetry of 
the surface. For the (110) surface we also evaluated the 
strain derivative of the surface stress by compressing the 
surface isotropically up to 12 %. The surface energy of 
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FIG. 5. The equilibrium crystal shape of InAs in (a) In-rich 
and (b) As-rich environment. The Miller indices of the sur- 
faces are noted. 

the strained surface is computed by subtracting the cor- 
responding energy of strained bulk from the total energy 
of the relaxed strained slab. The surface energy as a 
function of the strain is fitted to a second order polyno- 
mial. The second order coefficient gives the sum of the 
derivatives of the surface stress + 26'ii22 + *5'2222 

which is equal to —0.5 eV/A^. Straining the InAs (110) 
surface epitaxially to the GaAs lattice constant would re- 
sult in a contribution from the second order terms to the 
surface energy smaller than 2 meV/A^, while the linear 
term amounts to ~ 6 me V/A^. Thus we neglect in the 
following the second and higher order corrections to the 
surface energy. 

Under thermal equilibrium conditions, the shape of a 
large InAs crystallite is given by the condition of lowest 
free energy. This can be obtained for zero temperature 
by applying the Wulff construction using the surface en- 
ergies of Tab. |ll[ The resulting shape is shown in Fig. 
|. As this fi gure is constructed using only the surface 
energies of the {110}, {100}, {111}, and {111} orienta- 
tions, we note that in general also higher index surfaces 
might be present, but the low Miller- index surfaces are 
expected to remain clearly prominent. Fig. || shows that 
under As-rich conditions all four considered surface ori- 
entations coexist on the equilibrium crystal shape. This 
result is in agreement with the shape of large, and thus 
presumably fully relaxed, IpAs islands grown on a GaAs 
substrate by Steimetz et alSA using metal-organic vapor- 
phase epitaxy. Under In-rich condition the {111} facets 
do not exist because they are energetically unfavorable. 
This is probably due to the fact that we did not consider 
the (^19 X ^19) structure. 

IV. EQUILIBRIUM SHAPE OF ISOLATED 
QUANTUM DOTS 

The equilibrium shape of a strained coherent island of 
a given number of atoms is determined by the minimum 
of its total energy with respect to its shape. To deter- 
mine the optimum island shape as a function of volume 
we follow the procedure already outlined in Ref. ^ As 
described in Section I the total energy is evaluated by 




FIG. 6. The different island shapes which are investigated. 
They consist out of {101}, {001}, {HI}, and {111} facets. 



summing the elastic energy, the strain renormalized sur- 
face energy and the edge energy. 

Accurate values of the edge energies are not known. 
We examined the influence of edge energies by calculating 
the equilibrium island shapes assuming the same value of 
the edge energy for all types of edges. We found that the 
island shape was not influenced as long as edge energies 
are smaller than 100 meV/A for a quantum dot of 10,000 
atoms, mainly because the edge energies only scale with 
V^/^. Heller et alL^ measured the energies of steps on 
the GaAs (100) surface to be 4 meV/A and 13 meV/A for 
the two rdjfferent types of steps. Recently, Kratzer and 
ScheffleiES computed a value of 25 meV/A for the one 
type. Edge energies should be of comparable size or even 
smaller and thus do not play a role for the island shape. 
Therefore, we neglect the edge energies in the following 
analysis. 

Due to the different scaling properties of the elastic 
and surface energy, the optimum island shape depends on 
the volume. We consider all possible island shapes which 
have low-index surface facets. An overview is shown in 
Fig. |[ If surfaces with other orientations and a smaller 
slope than the {101} facets would appear on the islands, 
these surfaces would facet into the thermodynamically 
stable orientations {101}, {001}, {111}, and {111}. We 
calculate the total energies of islands bounded by {101}, 
{111}, and {111} facets. The filled symbols in Fig. ^ are 



5 




ESurface/y2/3 (meV/A^) 

FIG. 7. The elastic energy per volume ^^1^^*^^/^ versus 
the surface energy per area ^«urface^y2/3 InAs islands. 
The symbols refer to the shapes displayed in Fig. ^: Square: 
square based pyramid with four {101} facets. Diamond: 
square based pyramid with two {111} and two {111} facets. 
Triangles up: huts with two {111} and two {111} facets. 
Triangles down: square based {101} pyramids with {111} 
truncated edges. Dots: islands with four {101}, two {111}, 
and two {111} facets. The open circles denote the corre- 
sponding truncated islands which are connected by the full 
lines. The dashed line is the curve of constant total energy 
^elastic + ^surface that sclccts the equilibrium shape for the 
volume V = 2.14 x 10^ A^ 

obtained using results from full finite element calcula- 
tions, whereas the small dots denote the truncated mesa 
shaped islands where the elastic energies are derived from 
the simple analytical approximation. The results of the 
elastic calculations are combined with the ab initio sur- 
face energies. Because the side facets of the islands are 
strained we include the first order correction of the sur- 
face energy due to the strain. For this we use the ab 
initio stress tensors from Tab. || and combine it with the 
strain field at the surface from the finite element calcula- 
tion. We have used the surface energy of the unstrained 
InAs (100) surface for the almost fully relaxed (100) top 
facets of the quantum dots, and the surface energy of an 
isotropically strained wetting layer for the area covered 
by the quantum dots. The contribution of the elastic and 
strain renormalized surface energies are displayed in Fig. 
1^, both of them divided by the their respective scaling 
factors. 

In Fig. ^ the optimum island shape is determined by 
the point where the line of constant total energy touches 
the manifold of the island energies from below. There- 
fore, the equilibrium island shapes for all volumes are 
given by the lower envelope of the manifold of the is- 
land energies. In Fig. || the equilibrium island shapes 
determined by this method are shown for two different 
volumes. The elastic energy scales with the volume V 




FIG. 8. The equilibrium shape of a strained coherent InAs 
islands in As rich environment at two different volumes, (a) 
V ^ 2x 10^ (10,000 atoms), (h) V ^ 4 x 10^ A^ (20,000 
atoms) . 

whereas the strain renormalized surface energy increases 
with the volume like F^/^, because the strain field of 
the side facets is invariant to the scaling. Due to their 
scaling the surface energy dominates at small volume and 
the elastic energy gains more importance at large volume. 
Thus, larger islands are steeper than smaller ones. As the 
{111} and {111} facets are steeper than the {101} facets 
they become more prominent on larger islands. There- 
fore, the island shape is not fixed, but changes continu- 
ously with the volume. This change in the shapes implies 
that the simple scaling laws which are valid for a fixed 
shape do not apply. 

The influence of the surface stress is in such a way that 
surface energies of the {111} and {111} facets are further 
lowered as compared to the {101} facets. Roughly 30 % 
of the surface of each facets are strained while the rest is 
completely relaxed. This means that the {101} facets are 
lowered due to the strain by A7 = — 1.7 meV/A^ whereas 
the {111} and {111} facets are lowered by A7 = —2.1 
meV/A^ and A7 = —4.1 meV/A^, respectively, i.e., the 
{III} facets are lowered 6 % more in surface energy than 
the {101} facets. Thus, the {111} faces dominate even 
more due to the surface strain. Although influenced by 
strain effects the equilibrium island shape remains similar 
to the ECS. 

Recent atomic, force microscopy (AFM) studies by 
Georgsson et alu of uncapped (not overgrown) InP is- 
lands on a GalnP substrate corroborate our theoretical 
predictions. The experimental island shape displays the 
same facets as our calculated islaiid shape. Although InP 
differs from InAs, ab initio result£j yield a similar behav- 
ior for the surface energies of InP as for InAs. Therefore, 
we expect to obtain a similar theoretical shape for InP 
islands. For InAs, however, the experimentally observed 
shapes differ ^om ours and show various different shapes. 
Moison etpp/.B observe {410} to {110} facets and Ruvk 
mov et aln only {110} facets. Whereas Leonard et alu 
report their island shapes as planoconvex lenses with a 
radius to height aspect ratio of about two. The great di- 
versity of experimental island shapes and the difference 
to ours indicates that the respective growth conditions do 
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not represent thermodynamic equilibrium but are driven 
by kinetics. Thus, kinetic effects such as Jesson et a/.El 
suggested may play a role. These kinetic effects could not 
only affect the distribution of the island size but also the 
shape of the islands. However, it should be possible to 
achieve thermodynamic equilibrium by choosing appro- 
priate experimental conditions. Alloying of the quantum 
dots with the GaAs matrix will also pla^ a role and might 
affect the shape of the quantum dots.t^ 

To gain further insight into the shape of the islands, 
more experimental investigations should be performed 
such as high resolution STM or AFM of the side facets. 
This could help the theory for further investigations. 
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